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Abstract 

The dynamic structure factor of a trapped Bose-Einstein condensed gas is 
investigated at zero temperature in the framework of Bogoliubov theory. Dif- 
ferent values of momentum transfer are considered, ranging from the phonon 
to the single-particle regime. Various approximated schemes are discussed, 
including the local density approximation, where the system is locally de- 
scribed as a uniform gas, and the impulse approximation, where the response 
is fixed by the momentum distribution of the condensate. A comprehensive 
approach, based on the eikonal expansion, is presented. The predictions of 
theory are successfully compared with the results of recent two-photon Bragg 
scattering experiments, both at low and high momentum transfer. Relevant 
features of the dynamic structure factor are also discussed using the formal- 
ism of sum rules and the concept of scaling. Particular emphasis is given to 
the regime of high momentum transfer, where the dynamic structure factor 
is sensitive to the behaviour of the order parameter in momentum space, and 
some instructive examples showing the consequence of long-range coherence 
are presented. 

PACS numbers: 03.65.-w, 05.30.Jp, 32.80.-t, 67.40.-w 
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I. INTRODUCTION 



The dynamic structure factor provides an important characterization of the dynamic 
behaviour of quantum many-body systems. In particular, its exploration has played a crucial 
role in understanding the physics of superffuid 4 He, starting from the measurement of the 
roton spectrum [|lj until the more recent determinations of the condensate fraction available 
from neutron scattering experiments ||. The dynamic structure factor is measurable through 
inelastic scattering, in which the probe particle is weakly coupled to the many body system so 
that the scattering may be described within the Born approximation |J. In the case of dilute 
gases it can be measured via inelastic light scattering as recently shown by the experiments 
of Refs. HH carried out on a trapped Bose gas of sodium atoms. The dynamic structure 
factor provides information on both the spectrum of collective excitations which can be 
investigated at low momentum transfer, and the momentum distribution which characterizes 
the behaviour of the system at high momentum transfer, where the response is dominated 
by single-particle effects. 

In superffuid helium the typical momentum giving the transition between the collective 
and the single-particle behaviour is fixed by the inverse of the range of two-body forces, a 
value close to the average interatomic distance. At larger momenta one explores microscopic 
features of the system which are sensitive to short range correlations and to the details of 
the two-body interaction. The situation is very different in a dilute gas where the transition 
takes place at momenta much smaller than the inverse of the scattering length, which fixes 
the range of interactions. As a consequence, in a Bose gas one can explore a domain of 
relatively high momenta, where the response of the system is not affected by collective 
features, nor by short range correlations, but is determined by the momentum distribution 
of the condensate. While in a uniform system this distribution is a simple 5-function, in a 
trapped gas it exhibits a non trivial behaviour and is strongly affected by the presence of 
two-body interactions. 

The purpose of this paper is to provide a systematic theoretical discussion of the be- 
haviour of the dynamic structure factor of inhomogeneous Bose-Einstein condensates at 
zero temperature and to make quantitative comparisons with the recent data obtained with 
two photon Bragg scattering experiments @J^], pointing out the role of two-body interac- 
tions at both low and high momentum transfer. The applicability of both the local density 
(LDA) and of the impulse approximation (IA) will be discussed in detail. The LDA assumes 
that the system can be locally described as a uniform gas, and is adequate at moderately low 
values of momentum transfer. Conversely, the IA assumes that the dynamic structure factor 
is sensitive to the momentum distribution, which, for Bose-Einstein condensed systems, is 
determined not locally but globally, according to the size and shape of the condensate wave- 
function. A comprehensive description of both the LDA and IA regimes will be presented 
using the eikonal expansion, which holds in the single particle regime at high momentum 
transfer. Special emphasis will be given to sum rules as well as to the scaling behaviour 
exhibited by the dynamic structure factor in the IA regime. Finally we will point out the 
occurrence of interesting features exhibited by the dynamic structure factor in the presence 
of vortices and of interference effects in momentum space. 
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II. DYNAMIC STRUCTURE FACTOR AND BOGOLIUBOV THEORY 



The dynamic structure factor of a many body system is defined by the expression 

5(q, E) = 1 e ~ PEm I HP»\ 2 S(E -E m + E n ) (1) 

mn 

where q and E are the momentum and energy transferred by the probe to the sample. In Eq. 
(H) |n) and E n are the eigenstates and eigenvalues of the Hamiltonian of the system, e~^ En 
is the usual Boltzmann factor, p q = e iq ' r ^ is the Fourier transform of the one-body 
density operator, and Z is the usual canonical partition function. 

In Refs. MM, the dynamic structure factor of a trapped Bose-Einstein condensate is 
measured using two-photon optical Bragg spectroscopy. Two laser beams are impingent 
upon the condensate. The difference in the wavevectors of the beams defines the momentum 
transfer q, and the frequency difference between the beams defines the energy transfer E. 
The atoms exposed to these beams can undergo a stimulated light scattering event by 
absorbing a photon from one of the beams and emitting into the other. After exposure to 
these laser beams, the response of the condensate is measured by a time of flight technique 
by which the number of optically excited atoms can be counted. The momentum transfer in 
each of these experiments was fixed by the particular optical setup, while the energy transfer 
was scanned by varying the frequency difference between the beams. Because atoms could 
be scattered by absorbing a photon from either of the laser beams, the response of the system 
actually measures the combination S(q, E) — S{— q, —E), thus significantly suppressing the 
effects of temperature in the measured signal. 

Let us start our discussion by recalling that in an ideal uniform gas Eq.([l|) takes the 
simple integral form || 

S(q,E) = y , dp„( p )[i±n(p + q)]^-^^)! + ^. ( 2 ) 

where n(p) = (cJpOp) is the statistical average of the momentum operator aj,a p , and the 
sign + (— ) holds for Bose (Fermi) statistics. The scattering process is hence enhanced and 
suppressed in Bose and Fermi gases respectively. In the fermionic case this reflects the 
Pauli exclusion principle. For large momentum transfer q the exchange term of Eq. (0) is 
negligible because the momentum distribution decreases rapidly at high momenta and one 
finds the expression 

S„(4^ = /*'(s-^ + £)»(P>, (3) 

known as the impulse approximation (IA) jTj. It is important to stress that the validity of 
Eq.([3]) is not restricted to the ideal gas, but holds in general at high momentum transfer 
also for interacting and non uniform systems, independent of quantum statistics. Of course 
in this case the momentum distribution n(p) will differ significantly from that of the ideal 
gas. 

At small momentum transfer the ideal gas prediction (0) is inadequate, especially in 
the case of Bose gases, in which mean field interactions drastically modify the structure of 
5(q, E), giving rise to the propagation of phonons. 
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The calculation of S"(q, E) in interacting many-body systems requires in general a major 
theoretical effort. In the following we will limit ourselves to the T = case and to the 
study of dilute Bose gases where Bogoliubov theory is applicable. This restricts the range of 
momenta q to the "macroscopic" regime qa <C /i where a is the s-wave scattering length. For 
larger values of q short range correlations become important and Bogoliubov theory is no 
longer adequate. In the conditions of the experiment of [Q, carried out on a gas of sodium 
atoms, the Bogoliubov approach is well applicable since q/H ~ 20 /mi -1 and qa ~ 0.06 H. 
The macroscopic condition is even better satisfied in the experiment of where smaller 
values of q have been used. 

According to Bogoliubov theory the excited states of the system are given by the solution 
of the coupled equations f| 



eu r 



-ev r 



_^ V 2 + Kxt(r)-/i + 2^ 2 (r) 
2m 

" ^V 2 + Kxt(r)-/i + 2^ 2 (r) 
2m 



u i r ) + ^o( r M r ) 



(4) 
(5) 



for the "particle" and "hole" components of the elementary modes. In Eqs. (^-jD e is the 
energy of the excitation, and K x t(r) is the external potential for which, unless differently 
specified, we make the axially symmetric harmonic choice: 



V cxt (r) = ^mu 2 ± (x 2 + y 2 ) + ^muj 2 z z 2 . 



(6) 



Furthermore [i is the chemical potential, g = 4nH 2 a/m is the coupling constant which will be 
assumed to be positive, and finally ipo(r) is the order parameter characterizing the ground 
state of the system. 

In terms of the functions u n (r) and v n (r), which satisfy the ortho-normalization condition 



J dr[u* n {r)u m {r) - v*(r)v m (r)] = 5 nm , 
the relevant matrix element of the density operator takes the form 

(n|p q |0> = / dvK(v) + <(r)]e^ r /Vo(r) 



and the dynamic structure factor then becomes, at T = 0, [10 



S( q ,E) = J2 y"dr[<(p)+<( 



(r)]e^- r /Vo(r) 



5(E - e n ) . 



(7) 



(8) 



(9) 



In the case of a uniform gas the "particle" and "hole" components are plane waves: 
u(r) = U exp (ip ■ r/H) and v (r) = V"exp (ip • r/H) and the coupled Equations (|]-|) give rise 
to the Bogoliubov excitation spectrum fTT|j : 



e(p) 



/ p 2 ( p 2 
2m V 2m 



+ 2gn 



(10) 
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where n is the density of the gas. This spectrum exhibits a phonon dispersion e = cp at low 
momenta, with the velocity of sound given by c = y^gn/m, while in the opposite limit of 
high momenta it approaches the free particle energy e = p 2 /2m. The transition between the 
collective and the single particle behaviour occurs at momenta of the order of H/£, where 

y/Sirna 

is the so-called healing length. 

Using the normalization condition (|7|), one finds the result 

S B(^E) = N^- ) 5(E-e(q)) (12) 

for the dynamic structure factor, which consists of a ^-function centered at the Bogoliubov 
frequency ( |I0|) . Eq. fll2]) yields the Feynman-like result 

for the static structure factor which tends linearly to zero at low momenta, and saturates 
to 1 in the opposite, high q limit. 

Results dT^ - Tl3| ) hold for uniform Bose gases. In the presence of non uniform trapping a 



natural generalization is provided by the local density approximation (LDA), according to 
which the system behaves locally as a piece of uniform gas whose dynamic structure factor 
is given by the expression fll2] ), evaluated at the corresponding density |jT2| : 

S LDA (q, E) = fdv n{v)8{E - e(r, q)) - f . (14) 
J 2me(r, q) 

In Eq. ( |T4] ) n(r) is the ground state density of the system, and e(r, q) = e(n(r),q) is the 
local Bogoliubov dispersion (fLOf). Eq. ( |I~4"D is expected to describe accurately the dynamic 
structure factor for momenta larger than h/R where R is the radius of the condensate, since 
the effects of discretization in the excitation spectrum can be safely ignored. At the same 
time the momentum transfer q should not be too large since the local density approximation 
ignores the Doppler effect associated with the spreading of the momentum distribution of 
the condensate, which is expected to become the leading effect in the dynamic structure 
factor at very large values of q. The Doppler broadening is accounted for by the impulse 
approximation (Q), which, however, ignores the mean field effects of Bogoliubov theory. The 
conditions of applicability for both LDA and IA will be established in the next section. 



III. LOCAL DENSITY, IMPULSE AND EIKONAL APPROXIMATIONS 

In the previous section we introduced two useful approximations to the dynamic structure 
factor of an interacting Bose gas: the impulse (|3]) and local density (|14]) approximations. 
These two descriptions hold in different regimes of momentum transfer. The purpose of 
this section is to discuss the corresponding predictions and conditions of applicability. We 
will also present a comprehensive description of the high q response of the system, based on 
eikonal expansion, which includes the LDA and the IA as special cases. 
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A. Local Density Approximation 



Let us first discuss the local density approximation flUp. An explicit expression for 
5*(q, E) can be obtained working in the Thomas- Fermi limit Na/a^ Q ^> 1, where the ground 



state density is given by (see, for example, [|13|]) 



n(r) = -{u- Kxt(r)) 
9 



(15) 



and the chemical potential takes the form 



Na\ 2/ \ 
15 Tkj ho 



(16) 



In Eq. ( |T6"D a^o = a/ h/muj^o is the oscillator length calculated using the geometrical average 
uj ho = {tjj^ijJz) 1 ^ of the oscillator frequencies. Using this density profile and the result of 
Eq. (PUD for the excitation spectrum, one obtains S(q, E) for a trapped Bose condensate as 
ITT" 



S LDA (q,£) 



15 (E 2 - E 2 r \ 
8 E r fx 2 



(E 2 - E 2 ) 



1- 



2£ ry u 



(17) 



where 



2m 



(18) 



is the free recoil energy. Differently from the case of a uniform gas (see Eq.(|12|)), the dynamic 
structure factor is no longer a 5-function, its value being different from zero in the interval 
E T < E < E r \J\ + 2/j,/ ' E T . The value E = E T corresponds to the excitation energy in the 
region near the border where the gas is extremely dilute and hence non-interacting. The 
value E = E T ^1 + 2fi/E v is the excitation energy of a Bogoliubov gas evaluated at the 
central density. Notice that the LDA expression ([17]) for S*(q, E) does not depend on the 
direction of the vector q even in the presence of a deformed trap. 

As pointed out in the previous section, the local density approximation requires that 
the momentum transfer be larger than H/R. For smaller values of q the response of the 
system is sensitive to the discretized modes of the system and the LDA cannot be longer 
employed. The theoretical analysis of the dynamic structure factor in this regime of low 
momentum transfer has been carried out in ||10|| . In the following we will always assume 
that the condition q ^> H/R is satisfied. 

Starting from fllT]) one can evaluate the inelastic static structure factor S*(q) = 



N 1 J 5*(q, E) dE, which takes the analytic form [|J: 



w 4 1 4a 2 



(3 + 2a- a 2 
16a 5 / 2 



7r + 2 arctan 



a — 1 
2y/a~ 



(19) 



with a = 2{i/E r . When g — ^ Eq. ([OJ) vanishes as S(q) ~ q/(2mc), reflecting the role 
played by dynamic correlations which strongly suppress light scattering in the phonon regime 
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(see Fig. [I]). Here c = 32a/ /i/m/157r corresponds to an average sound velocity, which, as 
expected, is smaller than the value \J ' n/m calculated in the center of the trap. 

Useful information about the dynamic structure factor can be obtained by evaluating its 
energy moments m.fc(q) = J^E k S(q, E)dE. For example the average excitation energy E 
and the rms width A rms can be defined as 

B = n±, (20) 

m 

and 



' m 2 I irii ' 



A rms =W— - ^ (21) 
V m o V m o/ 

respectively. 

Let us first discuss the phonon regime in which E x <C /x, or, equivalently, g£ <C h, where 
£ is the healing length ( |TT| ) calculated at the central density. Typical values of £ in sodium 
samples, where a = 2.75 nm, are of 0.1 -j- 1 /im depending on the density of the gas. Using 
the LDA expression (IT7) for 5*(q, E) one finds, in the phonon regime, the results E = cq and 



A rms ~ 0.3 cq, showing that the width of the signal is not much smaller than the average 
energy. The shape of 5*(q, E) turns out to be asymmetric as a function of E (see Fig. [|), 
and the peak response occurs at an energy which is higher (by about 15%) than E. 

The dynamic structure factor S(q, E) has been recently measured in the phonon regime 
||. A typical experimental curve is reported in Fig. |2| together with the prediction (p"7j). 
These measurements also show (see Fig. H) the static structure factor S(q) in the phonon 
regime to be smaller than the non interacting gas value S(q) = 1, in agreement with the 
predictions of Eq. (0) @. 

In the opposite limit of large momentum transfer, where g£ ^> h, the excitation energy 
(p0|) predicted by the LDA is given by the expression 

E = E r + j fx . (22) 

This result corresponds to the average of the Bogoliubov energy e(q, r) = q 2 /2m + gn{v) 
holding at high q. Notice that, due to the asymmetric shape of the dynamic structure factor 
(|T7|), the average energy ( |2"2"D turns out to be smaller than the peak energy 

2 

-E^peak = E v + -/i . (23) 

This asymmetry should be, in principle, taken into account in the fit of experimental data. 
However, the shift of the line center was effectively determined from experiments in || using 
symmetric Gaussian fits and was in good agreement with Eq. ( [22] ) (see Fig. |3|). In the same 
regime of large q the rms width predicted by the LDA is given by 



A L » A = \l—ti (24) 



and is independent of q. 
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Thus, for momentum transfers q which are sufficiently small to allow the use of the local 
density approximation, the structure factor of an inhomogeneous Bose-Einstein condensate 
can be derived from the Bogoliubov spectrum for a uniform condensed gas. Studies in this 
regime therefore serve as a probe of the spectrum of both collective (q£ -C K) and free-particle 

3> h) excitations. 



B. Impulse Approximation 

Let us now discuss the response of a trapped condensate to very large momentum trans- 
fers at which the form of the dynamic structure factor is dominated by Doppler broadening. 
In this regime, the dynamic structure factor is correctly described by the impulse approx- 
imation (0). Inelastic scattering at such high momentum transfers allows one to directly 
measure the momentum distribution of a trapped Bose gas. The possibility of such mea- 
surements is highly appealing since most of experimental investigations in these systems have 
been so far limited to the study of density profiles. In current experiments on harmonically- 
confined Bose gases, the sizes of the condensate R and of the thermal cloud Rt are typically 
comparable. In the Thomas-Fermi regime, the ratio between the two radii is given as 

R I li /fkj ho f 15 N o a \ (25) 



Rt V k B T V k B T \ a ho 

where Nq is the number of atoms in the condensate. Due to the large value of the Thomas- 
Fermi parameter Noa/a^ Q , this ratio is typically close to unity. Expression ( p5| ) also provides 
an estimate for the ratio of sizes of the two components measured in time-of-flight exper- 
iments, in which the trap is suddenly switched off and the gas allowed to freely expand. 
While the expansion of the thermal cloud is indicative of the non-condensate momentum 
distribution before release from the trap, the expansion of the condensate in the Thomas- 
Fermi regime is dominated by the release of interaction energy and does not reveal its initial 
momentum distribution. 

In contrast, the distinction between the condensate and the thermal cloud in momentum 
space is stark. A confined condensate of finite size has a momentum distribution of width 
Ap c ~ H/R fixed by the inverse of the size R of the condensate. The momentum width of 
the thermal cloud is instead given by the temperature of the gas as Apx ~ \Jmk B T . For 
harmonic confinement in the Thomas- Fermi regime, one then finds 

APc ^ ( 15^] . (26) 



Ap T V k B T \ a ho 

In contrast with the comparison of the condensate and the thermal cloud in coordinate 
space (|25|) , the distinction between the two components in momentum space is strongly 
enhanced by two body interactions as the Thomas- Fermi parameter N a/a^ o increases. The 
investigation of the momentum distribution consequently provides a deeper understanding of 
the phenomenon of BEC. In particular the smallness of the width Ap c reflects the presence of 
long-range coherence. Measuring the momentum distribution at zero temperature would, in 
principle, give access also to the quantum depletion of the condensate. In practice, however, 



S 



the quantum depletion is too small and broadened over too large a momentum range [TJj to 
be observable in present experiments. 

The expression (H) for the impulse approximation can be also written in the form 



(27) 



where we have assumed that the vector q is oriented along the x axis, and p x = m(E — E r )/q. 
The integral fdp y dp z n(p x ,p y ,p z ) is also called the longitudinal momentum distribution. 
Eqs. (H) and (|27|) show that in the regime of applicability of the IA one can extract useful 
information on the momentum distribution starting from the experimental measurement of 
the dynamic structure factor. 

In a dilute Bose gas at zero temperature the momentum distribution is given by n(p) = 
| 2 , where 



[2nH) 



-3/2 



drif)(r) 



ip-r/H 



(28) 



is the Fourier transform of the order parameter. The form of n(p) for a trapped condensate 
has been discussed previously |l6,17j. In the Thomas- Fermi limit Na/a\ lQ ^> 1 one finds the 
simple analytic result 



(29) 



m 


IMP)] 




p2 



where J2(z) is the usual Bessel function of order 2, 

Na\ 1/5 



Ra 



15—1 

0-ho J 



A l/3 



(30) 



is 



the Thomas-Fermi radius of the condensate in the x — y plane, and p = 
\jp 2 x + Py + {Pz/^) 2 R±/fi is a dimensionless variable, with the parameter A = u z /uj± fix- 
ing the anisotropy of the external potential. 

Eq. (|29"D explicitly shows that the momentum distribution scales as l/R± and is conse- 
quently much narrower than that of the non-interacting gas 



wibg(p) = N 



exp 



a 



J^l/3 

n z 



pI+pI + j 



(31) 



since R± ^> a ho . 

In the impulse approximation (p7 ) the peak of ^(q, E) coincides with the recoil energy 
E T while the curve is broadened due to the Doppler effect in the momentum distribution. 
A useful estimate of the broadening can be obtained carrying out a Gaussian expansion in 
the dynamic structure factor ( p7|) near the peak value E = E T . One finds: 



Si A (q, E) ~ ^(q, E r ) exp 



{E-E r 
2A? A 



(32) 
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with Aj A = — [Si A (q, E) / d^S\pS,% £')]|e=e i - By calculating the second derivative of 
with the Thomas- Fermi profile ( p9l) for the momentum distribution, we obtain, after some 
straightforward algebra, the result 



A 



qH 
3 mR\ 



(33) 



The Gaussian profile ( |3~2D reproduces very well the exact curve (see Fig. ||), so that the 
Doppler width ( |33| ) can be usefully compared with experiments, where the widths are usually 
extracted through Gaussian fits to the measured signal. The Doppler width (^) is linear 
in q, and for large momentum transfer it can become comparable or even larger than the 
mean field width fl24|) . In Fig. ^ we show the theoretical prediction (|33|) together with the 
experimental values obtained at several densities. This figure confirms that the IA accounts 
for the observed widths in the low density regime. At higher density the mean field effect 
(H]) can no longer be neglected. 

It is worth noticing that the width (|33| ) should not be confused with the rms width ( |2T|) 
which requires the evaluation of the moment and, using the IA expression (§) or 
takes the form: 



A, 



rn 



-^kin > 



(34) 



instead of (|33|) . Here E£ in = J dp p 2 ,n(p)/2m is the x-component of the kinetic energy of the 



condensate. The evaluation of the kinetic energy requires a careful analysis [0,0 of the 
region near the boundary of the condensate, and cannot be evaluated using the Thomas- 
Fermi expression fl29f ) for n(p), which incorrectly yields a divergent result. For large N 
samples one finds [TM], assuming isotropic trapping, 



5h 2 



pa; 
^kin 



N 6mR 2 



log 



R 



1.3a ho 



(35) 



The logarithmic term reflects the fact that the m 2 moment, and hence the rms width, is 
sensitive to the high energy tails of the dynamic structure factor, a region which is difficult to 
measure since the intensity of the signal in the tails is very small. Because of this, estimate 
(|3~3|) is much more significant from the experimental point of view than expression (|34]). 

The investigation of the dynamic structure factor also provides information on the coher- 
ence effects exhibited by the system and in particular on the behaviour of the off-diagonal 
one-body density 



p( 1 )(s) = N J dRdv 2 ...dv N iP* (R+^,r 2 ,...,r N ^ i> (r - |, r 2 , r^) 

dpn(p) exp 



P ■ s 



(36) 



where ip(rx, rjv) is the many-body wave function of the system, and n(p) is the momentum 
distribution. By taking the Fourier transform of Eq. (E7|) with respect to p x , one finds the 
result 



p«(^,0,0) = JdES 1A (<i,E)exp 



-i^{E-E x ) 



hq 



(37) 
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which shows that the one-body density is a measurable quantity if one works at high q where 
S(q, E) ~ SiA(q, E). In a uniform Bose-Einstein gas p^(s x , 0, 0) tends to a constant when 
s x is large. In a finite system p^(s x , 0, 0) always tends to zero when s x — > oo. The typical 
length over which decreases can be of the order of the size of the sample or smaller 
depending on the degree of coherence. Using, for example, the Gaussian profile fl32|) for 
S(q, E) one finds 

p {1 \s x , 0,0) ~ iVexp 

with Xx = R± a/3/8. One can see from Eq. (^) that Xx plays the role of a coherence length 
IP (]|| , which turns out to be of the order of the size of the system. This result reflects the 
fact that in a Bose-Einstein condensate the Heisenberg inequality ARAp c > h/2 is close to 
an identity. 



(38) 



C. Eikonal Expansion 



In order to describe the transition between the LDA and IA regimes discussed above and 
to better understand the corresponding conditions of applicability, it is useful to evaluate the 
high energy solutions of the Bogoliubov equations (f§-|5|) using an eikonal expansion pl|,|22 



In the large q limit, where we are interested in the solutions with energy e much larger than 
the chemical potential /i, one can neglect the function v (r) in Eq.(||) and look for a solution 
of the form u(r) = exp [ipf ■ r/fo}u(r), where p/ is the momentum of the excitation and u(r) 
is a slowly varying function. Keeping only terms with first spatial derivative of u(r) (eikonal 
approximation), the solution of @ with energy e = p^/lm takes the form: 



u(r) ~ exp 



P/ 



exp 



rn 



—i- 



p f H 



dx'V e 



ea{x ,y,z) 



(39) 



where the effective potential V e s(r), calculated in the Thomas Fermi limit, is equal to gn(r) 
inside and to K xt (r) — p outside the condensate. At high q the main contribution to the 
dynamic structure factor (Q) arises from the excited states with p/ ~ q. This has been taken 
into account in the eikonal correction u(r) (second factor of Eq. (p9|) ), where p/ was chosen 
along the a;-axis, i.e. the axis fixed by the vector q. Notice that in the eikonal approximation 
the free particle solution e ipf ' r ^ n is modified by the interactions only through a change of 
the phase. The importance of such a correction in the behaviour of the dynamic structure 
factor depends on the maximum phase deviation of u(r) from a pure plane wave, which is 
determined by the Born parameter b 



p_qRj 

E r h 



(40) 



where E r = q 2 /2m and R is the Thomas- Fermi radius ([30]). Differently from the ratio p/E r , 
the Born parameter depends explicitly on the size of the atomic cloud. 
From Eq. (|5|) one gets, in first approximation, 



v r 



mgn[r) 
p) 



-u(r) 



(41) 
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Notice that inclusion of v(r) in Eq. (^) for u(r) would result in a higher order correction. 
Making the transformation P/ = p + q, the dynamic structure factor (^) takes the form: 

5(q, E) = [dp\M\*5 (e - $- - ^ ) , (42) 



2m m 



with 



M = (2^ S^ ^ ^ 



TfL 

dx'V P f 



o 



gn{r) 
q 2 /m 



exp 



P • r 

n 



(43) 

where we have approximated pf ~ g in the evaluation of v(r) and in the eikonal correction 
and neglected the kinetic energy term p 2 /2m in the argument of the 5-function. This is 
a very accurate approximation if one works with large samples, where the Thomas-Fermi 

approximation ^o( r ) = \J ~ Kxt( r )) 1 9 applies. 

If the Born parameter ( [40|) is small, then the eikonal correction can be neglected and, 
ignoring the small term gn(r)m/q 2 in Eq. (B3), one recovers the IA result (|3|). Conversely, 
if b is large one finds a different behaviour. In this case the main contribution to the 
double integral \M\ 2 arises from the region where \x\ — a^l ~ R±/b <C R± and the eikonal 
correction to the relative phase (m/Hq) f^dx'V e ff(x' , y, z) can be consequently written as 
(m/Hq)V e s(x, y, z){x\ — X2), with x — (xi + x-i)/2. In these expressions we have set y\ = 
y2 = y and z\ = Zi = z as a consequence of the integration on p y and p z in Eq. (^2|). 
By integrating with respect to the relative variable X\ — x<i and to p x one finally recovers 
the LDA result ( |T4] ) for large q, where e(q, r) ~ q 2 /2m + gn(r). In conclusion the eikonal 



approximation ( [42| - f43D provides the proper description of the dynamic structure factor in 
the Thomas-Fermi limit in all the regimes of momentum transfer ranging from the LDA to 
the IA, provided /i<c£ r . It is interesting to notice that the Born parameter (f40D fixes also 
the ratio between the widths (U) and (|33|) of the dynamic structure factor calculated in the 
IA and LDA limits respectively. In fact one has 

^ = - (44) 
A IA 14 ' ^ 

so that the comparison between the two widths provides an equivalent criterion for the 
applicability of the two opposite approximations. The transition between the LDA and the 
IA takes place when the ratio (|44]) is close to unity. Using the relation fi = h 2 /2mC, 2 , where 
£ is the healing length ( |TTD calculated at the center of the trap, this corresponds to the value 
q = HR±/1A^ 2 for the momentum transfer. Notice that this value is much larger than the 
inverse of £, since in the Thomas- Fermi regime R± 3> £ (see also Fig. 11). For example, 
in the cloud of sodium atoms explored in Ref. [ffl with peak density n(0) = 3.8 10 14 cm -3 , 
corresponding to b ~ 14, one has 1/R± — 0.1 /an -1 , l/£ ~ 5 /im _1 , and R±/1A^ 2 ~ 20 
jum -1 . 

Eqs. (|42] - |43|) can be easily calculated numerically in all the regimes between the LDA 
and the IA. The comparison with the available experimental results (see Fig. ^J) is rather 
good and explains the deviation of the observed signal from the LDA as well from the IA 
predictions. The width of the dynamic structure factor in general is well reproduced by the 
quadrature expression a/ ^Iba + ^ia accounting for both the LDA and the IA widths (see 
also pl). 
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IV. SUM RULES 



The conditions of applicability of the local density and impulse approximations presented 
in the previous section can also be discussed using a sum rule approach [|4|] which allows 
for an exact determination of the width of the dynamic structure factor in some relevant 
limiting cases. 

We have already introduced in the previous section the moments of the dynamic structure 
factor relative to the operator p q . In terms of the matrix elements (n|p q |0) of the density 
operator one can write 



m k { q) = y^\{n\p^)\\E n -E G ) k = dES(q,E)E k , (45) 

and, using the closure relationship ^2 n \n)(n\ = 1, one can easily express the moments 
m^(q) in terms of mean values on the ground state |0) and of commutators between the 
Hamiltonian and the operator p q . Using the property 5*(q, E) = S(—q,E), holding in the 
presence of parity or time-reversal invariance, we find the following results for the lowest 
moments: 



m (q) = (p{p (l )-\(p (i )\ 2 = NS(q) 



m 2 ( q ) = ([plH][H, Pci ])=N 
1 



m 3 (q) 



2 
N 



[2 - %)] 
([[plH},[H, [H, Pci ]]]) = 



2m 



h 2 a 2 



m c 



2m 



+ 4 



2m 



Tpx 



+ 



E i 



int 



N 



a 2 ft 2 
2m m 



(46) 
(47) 

(48) 
(49) 



where we took the vector q along the x axis. Notice that in Eq. fl46|) we have subtracted 
the elastic contribution |(p q )| 2 . The kinetic structure function D x (q) is defined by 



D x {q) = Jj J dridr 2 cos [q( Xl - x 2 )\ V?Vfp (2) (ri, r 2 ; r' 1; r' 2 )\ ri = r > it 



(50) 



where 

p (2) (ri, r 2 ; r[, r' 2 ) = N(N - 1) J dr 3 dr i ...dr N tp*(r 1 , r 2 , r A r)^(r' 1 , r' 2 , r 3 , r 4 , r N ) (51) 
is the two-body density matrix. In the asymptotic limit q — > oo this function is related to 



the kinetic energy of the system |25 



lim D*{q) = |^ 



(52) 



The 777,3 sum rule 
Hamiltonian 



has been obtained evaluating the commutators with the effective 



(53) 



i<j 
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and using the corresponding ground state of Gross-Pitaevskii theory. In particular E- mt = 
g J drn(r) 2 /2 corresponds to the expectation value of the two body interaction energy. 

The f-sum rule ( f4"T|) is model independent [EH] and is satisfied by both the LDA and IA 
as can be explicitly shown by integrating the corresponding dynamic structure factors. The 
other sum rules are instead correctly reproduced only in suitable ranges of momenta, which 
thereby provide the corresponding regimes of applicability of the two approximations. 

The non energy weighted sum rule fl4"E| ) coincides with the inelastic static structure factor 
S*(q), a quantity of high interest, directly related to the Fourier transform of the diagonal 
two-body density matrix: 

5(q) = 1 + ^ J dr.dvi e ^~^ [^(r,, r 2 ; r x , r 2 ) - n(n)n(r 2 )] . (54) 

The LDA prediction for 5(q) has been already discussed in the previous section [see Eq. (|T9|)] . 
This prediction is expected to hold in all regimes of small and high q except, of course, when 
qR < h. In fact the static structure factor 5*(q) is not sensitive to the Doppler broadening 
which instead affects other moments of the dynamic structure factor. An interesting property 
of the static structure factor is the occurrence of a 1/q 2 correction to the large q asymptotic 
value: 

= (55) 



7 q 2 

Such a correction is a peculiarity of dilute Bose gases [|Tj| and is directly related to the shift 
of the average excitation energy given by the Feynman ratio ni\ /tuq = E T + 4///7. The IA 
does not instead predict any g-dependence for the static form factor and consequently fails 
in reproducing the shift of the peak. 

The m2 sum rule is also interesting for understanding the difference between the LDA 
and the IA. For large q this sum rule contains two corrections to the leading asymptotic 
value: 



lim 7772(g) 

q~ >oo 



2m 



+ f 

2m 



r n 



(56) 



The first correction, fixed by the chemical potential, arises again from the large q behaviour 
(|55l) of S(q) which enters the expression (|48|) for m 2 . The second contribution, proportional 
to the kinetic energy, arises form the kinetic structure factor D x (q). The first correction is 
correctly given by the LDA, the latter by the IA. It is worth noticing that if one calculates the 
rms width (|2ll) only the kinetic energy term survives in the large q limit. This confirms the 
correctness of the impulse approximation in reproducing the width of the dynamic structure 
factor at high q. 

Finally the 7773 sum rule is interesting because it can be explicitly evaluated for any 
value of q. For a uniform gas (V exk = and E kin = 0) the ratio a/ m 3 /mi coincides with the 
Bogoliubov excitation spectrum (0). In the presence of harmonic trapping it is instructive to 
calculate Eq. (|49|) in the small momentum transfer limit, where only the last term, containing 
the external potential, survives and the ratio \[m?Jm\ coincides with the frequency uu± of 
the dipole mode. Indeed, the dipole mode is the only mode excited by the density operator 
Pq = J2j e iq -' T ^ H in the q — > limit. 
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Before concluding this section we stress that the results discussed above hold in the Bo- 
goliubov regime qa <C ft. For larger momenta, relevant for example in the case of superfluid 
helium, a different behaviour takes place. For example, result ([55]) for the static structure 
factor S(q) is no longer valid. Particular attention should be also paid to the kinetic energy 
which characterizes the large q behaviour of 777.2 and hence of the rms width. One should in 
fact distinguish between the kinetic energy of the condensate and the full kinetic energy of 
the system which, even in dilute Bose gases, is dominated at T = by the high momentum 
components of n(p). The former is given by (|35|) and becomes smaller and smaller as R 
increases. The latter is instead of the order of the interaction energy and would determine 
the rms width of the dynamic structure factor for momentum transfer larger than h/a. 



V. SCALING AND IMPULSE APPROXIMATION 

In the previous sections we have often referred to the impulse approximation as the proper 
theory to describe the high q response of the system. At the same time we have pointed out 
that the IA does not account for the mean field shift ( P^D of the peak energy occurring at 
high q, which is instead correctly predicted by the LDA. In this section we discuss in what 
sense the IA provides the exact asymptotic description of the dynamic structure factor. The 
discussion is simplified by using scaling, a concept already employed in other many body 
systems, including atomic nuclei [p7| , liquids and solids Let us introduce the scaling 

variable 

which is the relevant variable to describe the asymptotic behaviour of the dynamic structure 
factor. We define the scaling function F (Y) according to the asymptotic behaviour 

Fo{Y ) = Urn JL^M) , (58) 

where, in the limit, the excitation energy E varies with q in order to keep the value of 
Y fixed. Comparison with Eq. fl27|) shows that F (Y) coincides with the longitudinal 
momentum distribution 

F (Y) = 1 J dp y d Pz n{Y,p y ,p z ) • (59) 
In terms of the scaling function F (Y) the rms width takes the form 



\ ms = -J [°° dYF (Y)Y 2 . (60) 
m y J -00 

Furthermore one has J^dYFoiY) = 1, and J^dYY F (Y) = 0. 

The scaling result (|58| - |59D hold for a wide class of many body systems interacting with 
realistic two-body potentials, and is not restricted to dilute gases (see for example J^J and 
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references therein). This means that corrections to the IA, due to final state interactions, 
give a vanishing contribution to the limit fl58|) . In general the following expansion holds: 



qS(q,E) 



m 



m 



N 



F (Y) + -F 1 (Y) 
Q 



F 2 {Y) + 



(61) 



In the presence of parity or time reversal symmetry the scaling function F (Y) is symmetric: 
F (Y) = Fq(—Y). Conversely one has: Fi(Y) = —Fi(—Y). This suggests that the proper 
symmetrization [S(q, E) + S(q, 2E — E T )}/2 of the measured signal with respect to the recoil 
energy E T = q 2 /2m, would ensure a faster convergence to the scaling limit ||28|| . From an 
experimental point of view the direct verification of scaling, by changing q and E keeping 
Y fixed, is likely the safest criterion for checking the achievement of the IA regime. A 
peculiarity of dilute gases is that it is possible to reach the scaling regime for values of 
momenta where the Bogoliubov theory is still applicable. 

For a trapped Bose gas, the eikonal Eqs. ( f42| - f43"|) can be easily expanded for small values 
of the Born parameter (|40"D, corresponding to high values of q (b = 2mfiR±/foq). One finds: 



q S(q,E) 



777 



N 



R± 
H 



[fo(y) + bh(y) + ---] 



(62) 



where fo(y) and fi(y) are dimensionless functions of the variable y = (R±/H)Y, directly re- 
lated to the scaling functions of Eq. (|6l|): F (Y) = R±fo(y)/h and Fi(Y) = qR±bfi(y)/mH. 
The function fo(y) is given by: 



fo(y) = J^f d Py d Pz 



h^y 1 + Py + pj) 
y 2 + Pl+ P 2 Z 



(63) 



and is shown in Fig. [| together with the Gaussian expansion 



fo(y) 



15tt 
192 



exp 



— y 
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(64) 



yielding the results (|32| - |33D for the dynamic structure factor. The function fi(y) is given by 



h(y) 



1 7 \ 

dr 1 _r±J (r±p±) / dx\/l 
o Jo 



r\ sin {yx)x I 1 



x 



(65) 



and is also shown in Fig. |6|. 

Starting from Eq. (^) one can evaluate the shift of the peak with respect to the recoil 
energy E T due to the first correction to the IA. Imposing the condition c^S^q, E) = one 
finds, after some straightforward algebra, -E pea k = E T + 2/i/3, showing that result (^3|) for 
the line shift holds not only in the high q LDA regime where b is large, but also for small 
values of b. 
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VI. DYNAMIC STRUCTURE FACTOR AND VORTICES 



The study of vortices in trapped Bose gases is presently a challenging topics of both 
theoretical and experimental investigation. First experimental evidence of a vortices has 
been recently reported p0|j3l| . On the theoretical side the structure of vortices, the corre- 
sponding stability conditions, as well as their consequences on the dynamic behaviour of the 
condensate have already attracted the attention of many physicists. The identification of 
suitable methods of detection has also been the object of theoretical investigation. These 
include the expansion of the condensate [3^] , the shift of the collective excitation frequencies 
PB| and the occurrence of dislocations in the interference patterns |S3J . 

In this section we show that the measurement of the dynamic structure factor in the 
IA regime would represent a powerful tool to reveal vortices in a trapped Bose gas. In 
fact a vortex strongly affects the momentum distribution of the system. This can be easily 
understood by noting that the kinetic energy of a trapped condensate is roughly doubled 
by the addition of a vortex ||19|| . In contrast, the density distribution of a condensate in the 
Thomas-Fermi regime is only slightly modified by the vortex whose size is small compared 
to that of the gas. 

In the presence of a quantized vortex aligned along the z-axis the wave function of the 
condensate takes the form 



V>(r) = e llfi if} {r ± ,z) , 
where if)o(r±,z)) is the solution of the Gross-Pitaevskii equation PJ55] 
H 2 V 2 h 2 



m 



2m 2mr, 



+ -rr( UJ ± r ± + u z z ) + #o(r±, z) 



(66) 



(67) 



which contains the additional centrifugal term H 2 /2mr 2 _. Solutions of Eq. fl67|) have been 
obtained numerically in [ 36[| . The density distribution \ipo(r±, z)\ 2 exhibits a hole whose size 
is of the order of the healing length (pT|) of the gas £, which, in the Thomas-Fermi limit, 
is much smaller than the size of the condensate. Also in momentum space the distribution 
exhibits a hole as shown in Fig. [7]. This is the consequence of the phase in Eq. (|66|), which 
gives a vanishing value to the integral (p8|) at p± = 0, where p±_ is the radial component 
of the momentum vector p. The size of the hole is of the order of H/R, and consequently 
comparable to the total size of the condensate in momentum space. This can be easily seen 
calculating the momentum distribution in the Thomas- Fermi limit. In this limit the main 
effect of the vortex on the momentum distribution arises from the phase e lip , and one can 
safely use for ip the Thomas-Fermi expression ^/(/x — V cxt (r))/g, holding in the absence of 
the vortex. The result for n(p) can then be written in the form 



Wtf(p) = N 



R 3 ± 15 A 



2k 



d(p e 



lip 



dr±r± 



v 7 ^ 



dze ip \ 1 



where p is the scaled momentum vector p = (p x ,p y ,p z /X)R±/H, already introduced in sec. 
Ill B. Notice that in the Thomas-Fermi limit the effect of the vortex is factorized through a 
dimensionless integral. 
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In Fig. H we report the dynamic structure factor calculated in the IA (see Eq. (|3|)) 
with and without the vortex. The calculation was carried out for a gas of N = 10 4 87 Rb 
atoms trapped in a disk type geometry (A = a/8). For this low density sample the IA is 
very accurate. The double peak structure in S(q, E) reflects the occurrence of a peculiar 
Doppler effect, and represents a clear signature of the vortex. In fact the vortex generates 
a velocity field in the condensate with significant components both parallel and antiparallel 
to the momentum transfer q. 



VII. INTERFERENCE EFFECTS IN MOMENTUM SPACE 



Finally, let us discuss the dynamic structure factor in terms of the occurrence of interfer- 
ence phenomena in momentum space. Interference has been so far investigated in coordinate 



space by imaging two overlapping Bose Einstein condensates |37|]. However, even if the two 
condensates do not overlap in space, they can interfere in momentum space |38]. This opens 
the possibility of investigating interference while avoiding any interaction between the two 
condensates. 

Consider, for example, a double well potential and let us assume, for simplicity, that 
the condensates in the two wells (condensates a and b respectively) have the same number 
of atoms (see Fig. ||). If the distance d between the two wells is large enough to avoid 
overlapping, and if the potential acting on the condensates a and b can be obtained by a 
simple translation, then their wave functions can be written as: 



A( r ) = ( r + 



(69) 
(70) 



where ipo is the solution of the Gross Pitaevskii equation for each condensate. The Fourier 
transforms of Eqs. (|69| - |70D hence read: 



Up) = e^ d / 2 %( P ) 

h(p) = e-*» d / 2 Vo(p) 



(71) 
(72) 



having taken the displacement between the two wells along the x axis. Under the above 
conditions any linear combination 



Vc(r) =Va(r) + e^V6(r) 



(73) 



of the wavefunctions ( p9] - |70| ) corresponds to a solution of the Gross Pitaevskii equation. 
These combinations represent coherent configurations which exhibit interference patterns in 
the momentum distribution: 



n(p) = 2n (p) 



1 + cos ( + If 



(74) 
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where n (p) = |0o(p)| 2 - These patterns have interesting consequences on the shape of the 
dynamic structure factor which, in the IA, takes the form 



S(q,E) = 2S (q,E) 



1 



cos 



Yd 



where Y is the scaling variable (|57|) . The dynamic structure factor 
frequency period 



(75) 



5D exhibits fringes with 



Au 



AE 



q 

md 



(76) 



In Fig. [lC] we show a typical result for S(q, E) corresponding to a distance between the two 
condensates four times larger than their radial width. The position of the fringes depends 
crucially on the value of the relative phase between the two condensates. 



VIII. CONCLUSIONS 



In this paper we have provided a theoretical discussion of the dynamic structure factor 
>S(q, E) of a trapped Bose-Einstein condensate at low temperature. A first important aim 
was the development of the proper many-body formalism, based on Bogoliubov theory, to 
describe in a quantitative way the several interesting features exhibited by S"(q, E). These 
mainly concern the role of two-body interactions which sizably affect the response of the 
system in all the relevant regimes of momentum transfer. The possibility of providing 
accurate theoretical predictions for the dynamic structure factor is particularly appealing 
in view of the recent experimental data obtained via two-photon Bragg scattering. The 
available data are in general agreement with theory, and thereby provide a further important 
proof of the crucial role played by two-body interactions in these trapped Bose-Einstein 
condensed gases. Interactions affect the shape of S(q, E) at both small momentum transfer, 
where they are responsible for the propagation of phonons, and at high momentum where 
they show up in the shift of the peak with respect to the free recoil energy q 2 /2m as well as 
in the width which is sensitive to both mean field and Doppler effects. Various approximate 
schemes have been considered in order to better discuss the main physical features. These 
schemes are summarized in Fig. 11 where the importance of the characteristic length scales 
of the problem emerges clearly. At momentum transfer smaller than the inverse of the size of 
the system the response is characterized by the discretized normal modes of the system. This 
regime has been already discussed in other works and has not been considered here. At higher 
momenta the system behaves locally as a uniform gas. This is the range of applicability of the 
local density approximation (LDA) which successfully describes the excitation of phonons 
as well as, at momenta larger than the inverse of the healing length, the corrections to the 
free particle motion due to mean field interactions. At even higher momenta the response 
of the system cannot be longer described locally, because it is sensitive to the momentum 
distribution of the condensate, a quantity associated with long-range coherence effects. This 
regime is well described by the impulse approximation (IA), a theory currently employed 
to investigate the quasi-free response of various many-body systems. In our work we have 
limited the use of the IA to momenta smaller than the inverse of the scattering length, 
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i.e. to the range of applicability of Bogoliubov theory. In this regime only the momentum 
distribution of the condensate is relevant and one can safely ignore the higher momentum 
components of n(p) which would be crucial to describe the response at momentum transfers 
larger than H/a. An interesting feature emerging from our analysis is that the transition 
between the LDA and the IA regimes is characterized by an important physical parameter, 
the so-called Born parameter fl40|) which depends explicitly on the size of the system. The 
corresponding transition takes place at momenta fixed by the combination UR/^ 2 . We have 
shown that the transition can be accounted for by an eikonal treatment of the solutions of 
the Bogoliubov equations and a full calculation of S(q, E) has been presented in this regime, 
showing good agreement with experiments. Many of the relevant features exhibited by the 
dynamic structure factor in these trapped Bose gases have been also presented and discussed 
using the formalism of sum rules and the concept of scaling. 

A second important point emerging from our analysis was that the dynamic structure 
factor in the IA regime of high momentum transfer offers a new important investigation tool. 
Actually Bose-Einstein condensation in momentum space provides a deeper understanding 
of long-range coherence phenomena in comparison with the studies of the density profiles 
which have been so far the main object of investigation. In our paper we have discussed 
some significant examples of such opportunities, including the case of quantized vortices and 
of interference effects occurring in momentum space. In both cases we have predicted non- 
trivial features which should be visible in the dynamic structure factor as a consequence of 
the peculiar behaviour of the order parameter in momentum space. In the case of vortices the 
dynamic structure factor is characterized by a hole occurring at the recoil energy and whose 
size is comparable to the total width of the signal. In a second example we have shown 
that the dynamic structure factor of two spatially separated condensates should exhibit 
interference fringes which are the consequence of the coherence existing between the two 
condensates. 
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FIGURES 
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FIG. 1. Static structure factor S(q) at T = as a function of q£/h- The prediction ( |19|) of the 
local density approximation (solid line) is compared with the experimental points taken from ||. 
The result of the non-interacting model (dashed line) is also indicated. 
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perimental points are taken from Ref. ||. The trapping frequencies are u>± = 2tt 150 Hz and 
lo z = 2vr 18 Hz. 
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FIG. 3. Line shift as a function of the central density of the cloud. The straight line is the 

theoretical prediction (^), and the circles are the experimental points of Ref. 
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FIG. 4. Line width as a function of the central density of the sample. The solid line is the 
prediction of the eikonal approximation (42). The dashed and dotted lines correspond to the IA 
and LDA predictions (^) and (p4|) respectively. The theoretical results are compared with the 
experimental data of Q. Both the eikonal and the experimental values are obtained through the 
Gaussian fit S(q, E) exp [{E — E) 2 /2A 2 ] to the signal. The momentum transfer is taken along the 
x-axis. The trapping frequencies are lo± = 2ir 195 Hz and u) z = 2ir 17 Hz. 
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FIG. 5. Dynamic structure factor of a trapped Bose condensate at T = 0. The numerical 
predictions of the eikonal approximation (solid curve), IA (dashed curve) and LDA (dotted curve) 
are compared with the experimental data of Ref. normalized to reproduce S(q) = 1. The 
momentum transfer is taken along the x-axis. The trap parameters are the same as in Fig. 0. 
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FIG. 6. Dimensionless scaling functions in Thomas-Fermi regime: the fo(y) function (Eq. (|63|)) 
(solid line) is compared with its Gaussian expansion (^) (dotted line). The function fi(y) (Eq. 
(|65|)) is also shown (dashed line). 
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FIG. 7. Density profile (up) and momentum distribution (down) of a trapped Bose condensate 
as a function of the radial variables r±_ and p± , integrated along z and p z respectively, in the absence 
(dashed line) and in the presence of a quantized vortex (full line). These profiles correspond, in 
the absence of the vortex, to a central density n(0) ~ 1 10 14 cm -3 . The trap parameters are 
u) z = 2vr 220 Hz, and A = a/8. 
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FIG. 8. Dynamic structure factor of a rubidium trapped Bose condensate at T = in the 
presence (solid line) and in the absence (dashed line) of the vortex, calculated using the impulse 
approximation (|J). The momentum transfer is taken along the x-axis. The trap parameters are 
the same as in Fig. ^. 
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FIG. 9. Typical configuration to explore interference of two Bose condensates in momentum 
space: two harmonic wells are separated by a distance d = 20 [im; each well confines a condensate 
of sodium atoms with radial size R± = 5 /mi and central density n(0) ~ 0.7 10 14 cm -3 . 



0.045 



3 0.015 







1 1 

q — 21 hfim^ 1 \ 


i i 

\ 23 Na 
























1/ w . 











-15 -10 -5 5 10 15 

(E — E T )/h [kHz] 

FIG. 10. Dynamic structure factor relative to the configuration reported in Fig. ^ without 
(dashed line) and with (solid line) coherence between the two condensates (Eq. (f75|)). The mo- 
mentum transfer q is taken along the x-axis. The relative phase between the two condensates is 
chosen equal to zero. The trap parameters are the same as in Fig. ||. 
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FIG. 11. Relevant wave vector regimes characterizing the dynamic structure factor at T = 0. 
R is the size of the condensate (^) ; £ is the healing length ( |ll| ) . The ranges of applicability of the 
local density approximation (LDA) and of the impulse approximation (IA) are also schematically 
indicated. The momentum transfer corresponds to the condition Alda = Ata for the 

LDA and IA widths of the dynamic structure factor (see Eq. (|44|) ) . 
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